arXiv: 1505.03153vl [astro-ph.GA] 12 May 2015 


Astronomy & Astrophysics manuscript no. SFR_Z_final ©ESO2017 

January 25, 2017 


Central star formation and metallicity in CALIFA interacting 

galaxies 

J.K. Banera-Ballesteros 1 - 2 *, S. F. Sanchez 3 , B. Garcia-Lorenzo 1,2 , J. Falcon-Barroso 1 - 2 , D. Mast 4 , R. Garcia-Benito 5 , 
B. Husemann 6 , G. van de Ven 7 , J. Iglesias-Paramo 5,8 , F.F. Rosales-Ortega 9 , M.A. Perez-Torres 5 - 1() , I. Marquez 5 , 

C. Kehrig 5 , R. A. Marino 11 , J. M. Vilchez 5 , L. Galbany 12 ' 1 ’, A.R. Lopez-Sanchez 14, 15 , C. J. Walcher 16 , and 

the CALIFA collaboration 


(Affiliations can be found after the references) 
January 25, 2017 


ABSTRACT 

We use optical integral-field spectroscopic (IFS) data from 103 nearby galaxies at different stages of the merging event, from close 
pairs to merger remnants provided by the CALIFA survey, to study the impact of the interaction in the specific star formation and 
oxygen abundance on different galactic scales. To disentangle the effect of the interaction and merger from internal processes, we 
compared our results with a control sample of 80 non-interacting galaxies. We confirm the moderate enhancement (x 2-3 times) of 
specific star formation for interacting galaxies in central regions as reported by previous studies; however, the specific star formation 
is comparable when observed in extended regions. We find that control and interacting star-forming galaxies have similar oxygen 
abundances in their central regions, when normalized to their stellar masses. Oxygen abundances of these interacting galaxies seem to 
decrease compared to the control objects at the large aperture sizes measured in effective radius. Although the enhancement in central 
star formation and lower metallicities for interacting galaxies have been attributed to tidally induced inflows, our results suggest that 
other processes such as stellar feedback can contribute to the metal enrichment in interacting galaxies. 
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1. Introduction 

Interactions and mergers are identified as key mechanisms in in¬ 
creasing the star formation rate (SFR) in galaxies (e.g., Sanders 
& Mirabel 1996; Borne et al. 1999). In particular, luminous in¬ 
frared galaxies (LIRGs), which are associated almost exclusively 
to merger events, present strong episodes of star formation (e.g., 
Ellison et al. 2013b). The scenario proposed by these studies and 
hydrodynamic al numerical simulations (e.g., Mihos & Hernquist 
1996; Barnes & Hernquist 1996) suggests that this enhancement 
is the result of tidally induced inflows of gas that favors the in¬ 
crement of star formation activity in the central region of inter¬ 
acting and merging galaxies. On the other hand, large samples of 
interacting and merging galaxies in the Sloan Digital Sky Survey 
(SDSS, e.g., Ellison et al. 2008, 2013a; Scudder et al. 2012; Pat¬ 
ton et al. 2013) indicate that although central SFR is enhanced 
in these objects compared with non-interacting galaxies of sim¬ 
ilar stellar mass, on average this increment is rather moderate. 
This indicates that the triggering of starburst episodes depends 
on more parameters than the mere fact of galaxies being in inter¬ 
action (orbital configuration and intrinsic properties of the pro¬ 
genitors; e.g., Cox et al. 2006; Di Matteo et al. 2008). 

It is also worthwhile noting that statistical spectroscopic 
studies like the one above are performed by observing a fixed 
projected portion of the galaxies (e.g., SDSS or GAMA fiber 
sizes; Casteels et al. 2014; Robotham et al. 2014). This yields 
measurements on different galactic scales, making it difficult to 
assess the extension of star formation. In this regard, IFS ob¬ 
servations of interacting galaxies are crucial for understanding 
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whether the enhancement in the SFR is a global or a localized 
process. An example of the usefulness of such observations has 
been demonstrated by Wild et al. (2014). From a detailed study 
of the Mice galaxies, they found that there is a moderate incre¬ 
ment in the central SFR and no net enhancement in the global 
SFR. This study considers a single interacting system. Statis¬ 
tical studies of spatially resolved star formation in interacting 
galaxies are therefore required. Even more, to make a fair com¬ 
parison, a homogeneously observed sample of non-interacting 
galaxies is also required. Recently, numerical simulations are 
starting to explore the extend of star formation induced by the 
interaction. Moreno et al. (2015) find enhancement in the star 
formation on the central kpc scales and moderately suppressed 
activity at larger galacto-centric radii. 

In the above picture, these metal-poor gas inflows decrease 
the central metallicity in comparison to non-interacting galaxies 
(e.g., Kewley et al. 2006). As a consequence, merging galaxies 
are thought to contribute to the scatter in the mass-metallicity 
relation (hereafter M-Z relation) of star forming galaxies. How¬ 
ever, recent numerical simulations that included feedback pro¬ 
cesses have suggested that nuclear metallicity depends on more 
factors than the dynamics of the gas, such as the chemical en¬ 
richment due to the ongoing star formation, the stellar (and pos¬ 
sible AGN) feedback, and returned material by evolved stars into 
the interstellar medium (ISM) of the entire galaxy (Torrey et al. 
2012). The result of the interplay between these different pro¬ 
cesses in these simulations could lead, in some cases, to a de¬ 
pression in the central metallicity or even an enhancement with 
respect to isolated galaxies. Spatially resolved observational 
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studies are then required to shade some light on the chemical 
evolution of interacting or merging galaxies. 

Recently, IFS studies have been carried aimed at the under¬ 
standing of properties of the ionized gas in individual Hu re¬ 
gions for a large sample of galaxies. The radial gradient of oxy¬ 
gen abundance in galactic disks has been characterized, with an 
observed flattening in this gradient caused by interactions and 
mergers (Sanchez et al. 2014). They also confirmed the global 
and local nature of the M-Z relation (Sanchez et al. 2013). These 
studies have been possible thank to the integral-held spectro¬ 
scopic (IFS) data provided by the CALIFA survey (Sanchez et al. 
2012). CALIFA is an ongoing IFS survey aimed at obtaining 
spatially resolved physical properties, covering the full optical 
spatial extent (up to ~ 2.5 times the effective radius), of 600 
galaxies in the nearby universe (0.0005 < z < 0.03) of any mor¬ 
phological type. It is complete in a wide range of stellar masses 
( from ~ 10 9 5 to ~10 n M o , Walcher et al. 2014). This dataset 
allows us to study for the first time spatially resolved proper¬ 
ties of the ionized gas in a significant sample of interacting and 
merging galaxies to unveil the effect of interactions on different 
spatial scales and at different stages of the interaction. It also 
allows us to compare these properties with a sample of homoge¬ 
neous observations of non-interacting galaxies. 

In this article we are aiming to investigate the impact of inter¬ 
actions in the star formation and oxygen abundance on different 
galactic scales. To accomplish this, we study the emission-line 
flux maps extracted from the data cubes of interacting and post¬ 
merger galaxies. From these maps we derived integrated prop¬ 
erties, such as the Her equivalent width [EW(Ha) hereafter] and 
flux ratios essential for deriving oxygen abundances. This pa¬ 
per is organized as follows. In Section 2 we present the control 
and interacting samples, their ionized gas maps, and the aper¬ 
ture measurements. The comparison of these two samples on 
different galactic scales is presented in Section 3. Finally, we 
summarize our results in Section 4. 


2. Sample and data 

2.1. Interacting and control samples 

The galaxies that we use in this study are included in the CAL¬ 
IFA objects observed until May 2014. The selection of the in¬ 
teracting or merging sample is explained in a forthcoming arti¬ 
cle (Barrera-Ballesteros et al., submitted). Briefly, we first se¬ 
lected those objects with companions from the observed CAL¬ 
IFA sample that are within a projected distance of 160 kpc and 
systemic velocities lower than Av < 600 km s 1 . We noted that a 
large number of companion objects are not included in the CAL¬ 
IFA mother sample. To enhance the number of complete pairs 
in the survey, we conducted a complementary project aimed at 
observing companions not included in the original mother sam¬ 
ple (RI.s J.K. Barrera-Ballesteros and G. van de Ven). These 
observations and their reduction were carried out using the same 
strategy and pipeline as from the CALIFA survey. From the re¬ 
maining sample of galaxies without close companions, we visu¬ 
ally selected the objects with morphological evidence of merger 
event as post-merger systems. Such objects are morphologically 
disturbed, highly irregular, with evident tidal features or is there 
any signature of a recent merging event (Ellison et al. 2013a). 
The sample of interacting and merging consists of 103 objects. 
We compared our sample of interacting/merging galaxies with a 
control sample of 80 non-interacting galaxies to match the inter¬ 
acting sample in stellar mass and color. This sample is presented 


in Barrera-Ballesteros et al. (2014), along with a comparison of 
their stellar and ionized gas kinematics. 

2.2. CALIFA data cubes and emission-line extraction 

Galaxies presented in this study were observed at the Calar Alto 
observatory using the PPAK instrument (Roth et al. 2005). This 
instrument consists of 331 fibers of 277 diameter each, con¬ 
centrated in a single hegaxon bundle covering a field of view 
(FoV) of 74" x 64" with a filling factor of 60%. A three-point 
dithering technique is used to reach a filling factor of 100% 
across the entire FoV (see details in Husemann et al. 2013). For 
this study we used the data cubes with an intermediate nominal 
spectral resolution (V500 setup) of 4/A4 ~ 850 at ~5000A. The 
wavelength range of the V500 setup covers several bright emis¬ 
sion lines (3745-7300A): Ha, H [5, Hy, [On] 44 3728, 3726, 
[O hi] AA 4958, 5007, [Nn] AA 5755, 6584, [S ii] AA 6716, 6730, 
among others. The final data cube has an homogenized spectral 
resolution (FWHM) over the entire wavelength range of 6.0A, 
and the wavelength sampling per spaxel is 2.0A. The total expo¬ 
sure time per pointing is fixed for all the observed objects to 45 
min. The data reduction was performed by a pipeline designed 
specifically for this survey (see details in Sanchez et al. 2012; 
Husemann et al. 2013, and Garcfa-Benito et al. 2014 for an up¬ 
date to the second data release). 

To extract the flux from the ionized gas emission lines, we 
used the FIT3D package 1 (Sanchez et al. 2006, 2011). We ob¬ 
tain these fluxes for the following lines: Ha, H/3, [O m] AA 4959, 
5007, and [N n] 445755, 6584. The analysis, including errors, 
was performed on a spaxel-by-spaxel basis, which allows us to 
create flux maps for each of these emission lines, once subtracted 
the underlying stellar population (see an example of emission¬ 
line maps in Fig. 1). We have derived the EW(Ha) in each spec¬ 
trum by dividing the Ha flux by the average of the continuum 
emission adjacent to Ha emission line. This average continuum 
was derived by measuring the mean flux of the continuum in 
spectral windows of 30 A centered at blue-shifted and redshifted 
100 A away from the Ha emission line. 

2.3. EW(Ha) and oxygen abundance measurements 

Figure 1 shows an example of the wealth of the data provided by 
the data cubes. Even for a single galaxy, several structures can 
be observed in the emission line and the EW(Ha) maps, which 
is a proxy for the specific star formation rate (sSFR). However, 
in this study we limited ourselves to the integrated properties for 
each galaxy. Studies exploring the star formation and ionized 
gas two-dimensional distributions in interacting galaxies will be 
left for future works. Thus from the emission line flux maps, we 
obtained the integrated properties for each galaxy. We selected 
those spaxels with errors smaller than 10% of their flux within 
two apertures centered on the optical nucleus: a central and an 
extended aperture. 

To study the impact of the scale used to derive the properties 
of the ionized gas, we selected two different scales for the diam¬ 
eter of the central and extended apertures: a projected (5 and 30 
arcsecs) and the effective radius (0.3 and 2.5 R e ff). We derived 
the effective radius for each galaxy using the same methodol¬ 
ogy as outlined in the Appendix A of Sanchez et al. (2014). For 
the projected scale, we chose the size of the central aperture to 
be approximately twice the size of the spatial PSF in the sur- 

1 http://www.caha.es/sanchez/FIT3D/ 
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Fig. 1. Emission-lines properties of a galaxy in the control sample (top panels, NGC 0214) and in the interacting sample (bottom panels, 
UGC 312). Left panels show the Ha flux. White-solid and black-dashed circles represent the central and extended apertures in arcsecs and 
effective radius units, respectively. Middle panel shows the [Om],l 5007, emission line flux. The units in both flux maps are (loglO) 10 -16 erg s -1 
cm -2 arcsec -2 . Right panels show the logarithm of the absolute values of EW(Hor) in units of A. 


vey (FWHM ~ 2.5 arcsec). The central region in both samples 
covers between ~ 0.6 and ~ 3 kpc. The extended aperture was 
chosen to cover a wide portion of the galaxies. In terms of the 
Petrosian radius we covered a generous portion of the galaxies 
(i.e, 1.5 < rp < 8). The central aperture in units of R e ff is, in 
most of the cases, larger than the spatial PSF. It is also compa¬ 
rable to the central projected aperture. In most of the galaxies, 
the extended aperture in units of R e ff covers the entire FoV (see 
Fig. 1). In a few cases where the two companions in a binary 
system are covered by a single FoV, this aperture includes both 
galaxies (e.g., NGC 169). Our goal has been to compare how the 
properties of the ionized gas in the central and extended regions 
of interacting galaxies change when we compare it with a con¬ 
trol sample (see Secs. 3.1, 3.2, and 3.3). Nevertheless, our data 
also allow us to study the change of the ionized gas metallicity 
on different galactic scales (see Sec. 3.4). 

To obtain the EW(Ha) for each galaxy, we averaged the in¬ 
dividual spaxel absolute values included within each of the aper¬ 
tures. In the case of the observed oxygen abundance in each 
aperture, we used the empirical calibration 03N2 derived in 
Marino et al. (2013), which makes use of the logarithm of the 
fluxes ratios log ([Om] T5007/H/3) and log ([Nil] T6584/Ha). 
The 03N2 empirical calibrator has been proven to be a robust es¬ 
timator of the metallicity when it is no possible to have its direct 
estimation by auroral lines (Lopez-Sanchez & Esteban 2010). In 
Sec. 3.5 we study the impact of this oxygen abundance calibrator 
by contrasting our measurements using the above calibrator with 
the one presented by Pettini & Pagel (2004). 


3. Results 

3.1. BPT diagram 

We plot in Fig. 2 for each aperture and scale described in 
Sec. 2.3 the classical diagnostic diagram using the line ratios 
log ([O hi] 45007/H/?) vs log ([Nil] T5755/Ha) (Baldwin et al. 
1981, BPT diagram hereafter) in order to determine the frac¬ 
tion of star-forming galaxies in each of the samples for a given 
aperture size. We also included the demarcation line described 
by Kewley et al. (2001). Star-forming objects lie below this 
line, while AGN-powered sources are mostly located above it. 
On both size scales, the trends in each aperture are similar. For 
the central apertures (see left panels of Fig. 2) in both samples, 
a large number of galaxies are spread over the righthand side 
of the BPT diagram, from the star-forming to the AGN zone. 
However, some interacting objects lie in the so-called left branch 
(i.e., log ([N ii]/16583)/Ha) < - 0.6), associated with star-forming 
galaxies. In the extended apertures (right panels of Fig. 2), a 
larger percentage of galaxies lie in the star-forming zone, indi¬ 
cating that the integrated flux includes a larger number of Hu 
regions than in the central aperture. 

We note that for both aperture sizes in these BPT diagrams, 
there is a non-negligible number of interacting galaxies located 
in the left branch, whereas none of the control sample galaxies 
are located in this zone. As explained in Sanchez et al. (2014), 
objects located in this branch own a high percentage of young 
stellar population. For the central aperture, these galaxies are in¬ 
cluded in pairs or interacting systems (10 objects on both aper¬ 
ture scales). This also holds for the extended aperture (13 and 
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12 objects for the 30 arcsec and 2.5 effective radius apertures, 
respectively). Since these galaxies are expected to have a signif¬ 
icant percentage of young stellar population, these results sug¬ 
gest an increment in the SFR in interacting galaxies with respect 
to the control sample, particularly in binary systems. In the fol¬ 
lowing section, we study the impact on the change in the sSFR 
by the interaction/merger across the galaxies by means of the 
spatially resolved EW(Ha). 
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Fig. 2. BPT diagnostic diagrams for the control (blue points) and inter¬ 
acting (red points) samples measured within a central (left panels) and 
extended (right panels) apertures centered in the optical nucleus. Fluxes 
in top panels are derived using apertures on an arcsec scale, while fluxes 
in the bottom panel are obtained using an effective radius scale. The 
solid lines in each panel represent the division between star-forming 
and non-star-forming galaxies presented in Kewley et al. (2001). The 
typical uncertainty is plotted at the top of each panel. 


3.2. Central and extended EW(Ha) 

Once we distinguished between star and non-star forming galax¬ 
ies for each aperture size in both the control and interacting sam¬ 
ples, we explore how the EW(Ha) changes on the different scales 
of these galaxies. The subsample of star-forming galaxies is se¬ 
lected in both samples following the criteria used by Sanchez 
et al. (2014): objects with line ratios below the Kewley demar¬ 
cation line in the BPT diagnostic diagram in Fig. 2 and EW(Hcr) 
larger than 6 A. Our spatially resolved data allowed us to study 
the distribution of the EW(Hcr) at different aperture sizes in 
galaxies classified as star-forming and non-star-forming galaxies 
depending on the ionized gas fluxes in each aperture. From Fig. 2 
we note that non-star-forming sample includes mostly AGN-like 
objects, suggesting that the fraction of objects with a LINER-like 
spectrum is negligible. 

In Fig. 3 we plot the log(EW(Ha)) distributions of star¬ 
forming and non-star-forming galaxies for both samples at the 5 
and 30 arcsec apertures (EW centra i and EW exten d e d. respectively) 
for galaxies classified as star-forming and non-star-forming 
galaxies in the 5 arcsec aperture (see Sec. 2). In Fig. 4 we plot 
the same distributions using the integrated fluxes within the 30 


Central aperture selection (5 arcsec) 
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Fig. 3. Distributions of the integrated EW(Hff) for the interacting (red) 
and control (blue) samples using an aperture size of 5 arcsec. Top and 
bottom panels show the distributions of the star and non-star-forming 
galaxies, respectively. Left and right panels show the EW(Htt) dis¬ 
tributions measured using the 5 (EW cen trai) and 30 (EW ext ended) arcsec 
apertures, respectively. For this classification we used the integrated 
emission line fluxes ratios presented in Sec. 2 using the same selecting 
aperture size. 


arcsec aperture as classification. We repeat the same exercise 
using the 0.3 and 2.5 R cl[ aperture sizes, (see Figs. 5 and 6). 
We note that the EW(Ha) distributions in both samples using 
either the arcsec or effective radii scales are similar for the cen¬ 
tral and extended aperture sizes, respectively. To make a reason¬ 
able comparison with previous results using single-fiber spectro¬ 
scopic, we focus our analysis in the EW(H a) derived from the 
arcsec apertures. 

Although the fraction of galaxies selected in the central aper¬ 
ture as star forming is similar for the control and interacting sam¬ 
ples (33 and 35 objects, respectively), when we compare their 
EW centra i distributions, we find that a significant fraction of in¬ 
teracting galaxies (10/35) present higher values than does the 
entire control sample (see top left panel of Fig. 3). Most of these 
objects are either in close pairs or in binary systems with evident 
signatures of interaction (33/35). A majority of these interacting 
galaxies are late types (28/35). We also find clear statistical dif¬ 
ferences between these two samples. The median log(EW cen trai) 
of the star-forming interacting subsample is larger than the me¬ 
dian of the control sample ([1.41 + 0.07] log(A) and [1.09 + 
0.03] log(A), respectively) 2 . These interacting galaxies cover 
a wide range of logEW(Ho') in comparison to the control sam¬ 
ple (standard deviations of [0.38 + 0.07] log(A) and [0.39 + 
0.07] log(A), respectively). Moreover, a Kolmogorov-Smirnoff 
test (KS-test, hereafter) reveals that these two samples are not 
likely to be drawn from the same parent sample (p KS = 0.001). 
We transform the EW(Ho') to sSFR by means of the empirical 

2 Errors in the medians and standard deviations are obtained from a 
bootstrapping method. 
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Fig. 4. Distributions of the integrated EW(Hff) using an aperture size 
of 30 arcsec. As in Fig. 3 blue and red histograms represent control and 
interacting EW(Ho-) distributions, respectively. Top and bottom panels 
show the distributions for star-forming and non-star-forming galaxies, 
respectively. 


relation presented by Sanchez et al. (2013). For the interacting 
star-forming galaxies, the median sSFR is (2.3 + 0.2) x 10 10 
yr 1 , while the median sSFR for the control subsample is ap¬ 
proximately half the interacting value ([9.3 ± 0.9] X 10 11 yr -1 ). 

The distributions of EW extem jed for galaxies selected as star 
forming in the central aperture are fairly similar (see top right 
panel of Fig. 3), and median values are similar in both samples 
([1.07 ± 0.01] log(A) and [1.18 ± 0.02] log(A), respectively). 
The moderate sSFR enhancement we observe in central regions 
of interacting galaxies (~ 2.5 times) has been reported by pre¬ 
vious statistical studies (e.g., Kewley et al. 2006; Ellison et al. 
2008; Knapen & James 2009; Patton et al. 2013). However, we 
point out that those studies covered a portion that was larger than 
the central region. In fact, their aperture sizes are similar to our 
extended aperture where we find a moderate reduction of the 
sSFR compared to isolated galaxies (~ 0.74). 

In comparison to the star-forming galaxies, the non-star- 
forming galaxies selected by the central aperture in both samples 
present a rather similar distribution toward low EW centra i values 
(see bottom left panel of Fig. 3). We would like to point out here 
that these subsamples include objects located in the AGN zone 
in the BPT diagram or those with EW(Hor) smaller than 6 A. Al¬ 
though it is beyond the scope of this study, the similarity between 
these distributions deserve more attention so we can understand 
the impact of the interaction in triggering the nuclear activity of 
galaxies. We find however, a clear difference between the distri¬ 
butions of EWextended for the same sample of galaxies (see bottom 
right panel of Fig. 3). The median EW extenc j et i of the interacting 
sample is four times smaller than the control one ([0.20 + 0.08] 
log(A) and [0.82 + 0.05] log(A), respectively). 

The control star-forming galaxies selected using the ex¬ 
tended aperture almost doubles the interacting ones (see top left 


Fig. 5. Distributions of integrated EW(Ha) for the interacting (red) 
and control (blue) galaxies. The distribution of the panels is similar to 
Fig. 3. For these histograms we use an aperture size of 0.3 R e ff. 


Extended aperture selection (2.5 R eff ) 
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Fig. 6. Distributions of integrated EW(Ha) for the interacting (red) 
and control (blue) galaxies using an aperture size of 2.5 R c ff. The distri¬ 
bution of the panels is similar to Fig. 5. 


panel of Fig. 4). The median log(EW cen trai) for the interacting 
subsample is almost three times larger than the median of the 
control sample ([1.41 ± 0.1]log(A) and [0.99 + 0.03] log(A), 
respectively). For the same subsample of star-forming galax¬ 
ies, EWextended is similar for the interacting and control samples 
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([1.11 ±0.02] log(A), see top right panel in Fig. 4). AKS-test re¬ 
veals that these two samples could be drawn from the same par¬ 
ent sample {pks = 0.3). This suggests that the net enhancement 
of the sSFR occurs mainly in the central region of the interacting 
galaxies rather than in the outer regions. 

In the extended aperture, the selection of non-star-forming 
galaxies yields almost three times more interacting objects than 
control ones. The distribution, as well as the median of EW cen trah 
is similar for both samples ([0.18 ± 0.07] log(A), see bottom left 
panel in Fig. 4). On the other hand, EW ex tended present signifi¬ 
cant differences between the interacting and control samples (see 
bottom right panel in Fig. 4). The median EW ex t e nded for this in¬ 
teracting subsample is three times smaller than the control one 
([0.17 ± 0.06] log(A) and [0.64 ± 0.03] log(A), respectively). In 
fact, a KS-test indicates that these two distributions are not likely 
to be drawn from the same parent sample ( p KS = 0.001). 

Our analysis shows that regardless of the criteria we use to 
select star-forming galaxies, the interacting sample has a con¬ 
sistent moderate enhancement in the central sSFR when com¬ 
pared with isolated star-forming objects. On extended scales, 
both samples present similar distributions of EW exte nded and, de¬ 
pending on the aperture used to classify the star-forming galax¬ 
ies, similar or moderately suppressed total sSFR. Most of these 
interacting objects are late-type galaxies in pairs or merging sys¬ 
tems. These results suggest that moderate enhancement (~ 2-3 
times larger) in the sSFR only occurs in the central region of in¬ 
teracting galaxies. This scenario is consistent with gas been fun- 
neled to the central region of galaxies owing to the interaction: 
large amounts of gas (subsequently enhanced sSFR) is found in 
the central region of interacting galaxies, whereas the integrated 
sSFR remains similar or moderately suppressed. Our spatially 
resolved study agrees with the numerical simulations presented 
by Moreno et al. (2015). They find that the increment of star for¬ 
mation in merging galaxies is observed within the central kpc, 
while in outer galacto-centric radii, the activity is moderately 
suppressed. Our results also agree with studies of supernovae 
radial distributions (e.g., Habergham et al. 2012; Herrero-Illana 
et al. 2012) that find Type-Ibc supernovae, those more correlated 
to the ongoing star formation (Galbany et al. 2014), are more 
centrally concentrated in disturbed/interacting galaxies. 

In galaxies considered as non-star-forming, the central 
EW(Hor) distribution is similar between interacting and isolated 
objects, whereas the extended EW(Ho') is systematically smaller 
for the interacting sample. This points out that the process re¬ 
sponsible for the observed central EW(Ha) in non-star-forming 
galaxies for both samples could be similar. From the BPT di¬ 
agrams (see Fig. 2), this process could involve only AGN ac¬ 
tivity. On the other hand, low values of the extended EW(Ho') 
for interacting galaxies could imply either a low budget of gas 
available to be ionized or a weak radiation field across the merg¬ 
ing galaxies. Since these galaxies are the complement to the 
star-forming subsample, they cover a wide range of interaction 
stages and morphological types. A further division of the non- 
star-forming galaxies EW(Har) in the above properties led a few 
objects in each bin, making it difficult to derive reliable results. 
We require a large sample of merging galaxies in particular in the 
post-merger and remnant stages to derive statistical meaningful 
conclusions about the mechanisms responsible for the observed 
trends in the central and extended EW(Hcr) in non-star-forming 
merging galaxies. 
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Fig. 7. Left: Distributions of AO/H in the central (5 arcsec) aperture 
for the interacting (red) and the control (blue) star-forming galaxies, see 
text for details. Right: Distributions of AO/H using a central aperture 
of 0.3 R eff . 

3.3. Central oxygen abundances 

Observational studies suggest that the enhancement in the star 
formation activity in the central region of interacting galaxies is 
connected with a dilution of the oxygen abundance (e.g., Kewley 
et al. 2006; Ellison et al. 2008; Sanchez Almeida et al. 2014a, b), 
which in turns contributes to the scatter of the M-Z relation. This 
decrement is associated with the supply of metal-poor gas to 
the central region via inflows. To test this scenario, we stud¬ 
ied the scatter between the observed metallicity and the M-Z re¬ 
lation for a given stellar mass in star-forming galaxies. From 
the emission line fluxes, we derived the oxygen abundance of 
each galaxy classified in Sec. 3.2 as star-forming in the central 
aperture for both the interacting and control samples. Then, we 
provided the difference (or scatter) between the abundance de¬ 
termined from the flux ratios and the one expected from the M-Z 
relation (AO/H) derived by Sanchez et al. (2013). For the ob¬ 
jects included in the CALIFA survey, we used the stellar masses 
presented by Walcher et al. (2014). For the companions not 
included in the survey, stellar masses were derived in a simi¬ 
lar fashion to the CALIFA galaxies. As explained in (Sanchez 
et al. 2013), this empirical M-Z relation is derived for metal- 
licities at one effective radius. Nevertheless, we can transform 
this abundance at different radii, in particular to the central aper¬ 
ture radius, by means of the metallicity gradient in disk galaxies 
(Sanchez et al. 2014). 

In Fig. 7 we plot the distributions of AO/H for both the inter¬ 
acting and control samples in the central aperture (5 arcsec and 
0.3 R e ff, left and right histograms, respectively). The distribu¬ 
tion for the control sample in the 5 arcsec aperture is centered 
well on zero, with a median value of (0.001 +0.006) dex and a 
rather narrow distribution (<xAO/H control = 0.056 ± 0.006 dex). 
On the other hand, the median of the distribution from the in¬ 
teracting sample is slightly larger than the one from the control 
sample (0.03 + 0.01 dex). However, the interacting sample has 
a wider range of values in comparison with the control sample 
(crAO/H merger = 0.10 + 0.01 dex). A KS-test indicates that these 
two samples are likely to be drawn from the same parent sam¬ 
ple distribution ( p KS = 0.32). We find a similar distribution for 
the 0.3 R e ff aperture with median values of (-0.001 ± 0.006) dex 
and (0.021 ±0.02) dex for the control and interacting samples, 
respectively. The interacting sample is slightly more metal rich 
than the control sample (0.02 ± 0.02 dex). This indicates that in 
central regions, star-forming merging galaxies have similar oxy¬ 
gen abundance to the control sample. Even though both aper¬ 
tures provide similar results, the most reliable one comes from 
the 0.3 R e ff aperture size since it covers a similar physical region 
in all the galaxies in both samples. This first IFS census sug- 
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gests that even though there is an enhancement of the sSFR in 
the central region of interacting galaxies, the oxygen abundance 
remains similar to the one found in isolated galaxies. In other 
words, the scatter observed in the M-Z relation does not seem to 
be linked to the star formation rate of the galaxies presented in 
this study. 

Our results differ from the single-fiber spectroscopic stud¬ 
ies carried out with the SDSS galaxy-pair sample (Ellison et al. 
2008). They found a lower oxygen abundance of ~ 0.05 -0.1 dex 
between their galaxy-pair samples and control sample. One evi¬ 
dent reason for the disparity between these two studies could be 
attributed to the number of objects in each sample. SDSS pairs 
and control samples included 1716 galaxies in binary systems 
and 40095 galaxies, respectively, whereas our samples of star¬ 
forming galaxies include 33 and 35 objects in both interacting 
and control samples, respectively. Despite this difference in the 
size of the samples, we are able to account for the enhancement 
in the central SFR of interacting galaxies (see Figs. 3 and 5), in 
agreement with the SDSS studies, but not for the dilution in cen¬ 
tral metallicities. We note that our reported small variations in 
the metallicity in the central aperture may be biased by aperture 
effects or even the calibration used to derived the metallicity. In 
the following two sections, we explore the possible bias in our 
results for these two factors. 


3.4. Metallicity as function of the aperture size 

Our IFS data allows exploring the effect of measuring the gas 
metallicity for different aperture sizes and scales. In particular, 
we can study whether the lower abundances presented in single¬ 
fiber studies for the star-forming interacting galaxies (e.g., Elli¬ 
son et al. 2008) could be caused by aperture effects. More impor¬ 
tant, we will have an estimation of how the metallicity changes 
at different galactic scales. In Fig. 8 we plot the differences in 
the median metallicity value for the interacting and control star¬ 
forming galaxies at different aperture sizes. We plot these differ¬ 
ences for two aperture scales in arcsecs and in effective radius. 
The sample of galaxies selected for these plots are those classi¬ 
fied as star-forming galaxies in the central aperture of each scale. 

For the aperture scale in arcsecs, we observe that the differ¬ 
ence in metallicities between the two samples decreases as the 
aperture size increases, reaching a plateau at ~ 15 arcsec. Af¬ 
ter this size the differences in metallicity are rather constant for 
different aperture sizes (—0.01 dex). Using the effective radius 
of each galaxy as aperture scale, the difference in metallicity for 
the two samples decreases as the size of the aperture increases. 
For the largest aperture (2.5 R e fr), the interacting galaxies shows 
a lower metallicity of ~ 0.02 dex with respect to control galax¬ 
ies. Although there is a decrement in the metallicity of merging 
galaxies for large apertures, we do not observe the difference 
presented by previous spectroscopic studies (e.g., Kewley et al. 
2006; Ellison et al. 2008). Even more in the physical motivated 
aperture (R e ff), we find rather similar metallicities in the central 
regions with hints of a dilution of metallicities of the interacting 
galaxies on larger scales. This suggests that despite the inflows 
of ionized gas induced by the interaction, the metallicity in the 
central regions of interacting galaxies presents similar properties 
to those for isolated galaxies. On large scales, there seems to be 
a dilution in the metallicity of interacting objects. 


3.5. Metallicity as function of the abundance calibrator 

In general, the variations in metallicity between the control and 
interacting galaxies are not larger than 0.1 dex. This regime of 
small variations can be biased by several factors. In this sec¬ 
tion we explore and quantify the impact of these differences by 
the abundance calibrator. In Fig. 9 we plot the differences as a 
function of aperture sizes using the improved empirical calibra¬ 
tion for the 03N2 indicator given by Marino et al. (2013). We 
overplot these differences using a different calibration provided 
by Pettini & Pagel (2004). As we note above, the differences 
between the interacting and control sample are very small when 
using the improved calibrator, flattening out at ~ 15 arcsecs to 
—0.01 dex. However, when we use the calibrator from Pettini 
& Pagel (2004), we find that in larger apertures, the interact¬ 
ing galaxies have lower metallicities than control galaxies. Even 
more at ~ 20 arcsec the differences in metallicity reach the value 
found by Ellison et al. (2008). 

To explain this apparent difference between our results and 
those using another calibrator, we note that our sample is closer 
than the sample presented by Ellison et al. (2008). At the me¬ 
dian redshift of the SDSS-pair sample (z ~ 0.06, see their Fig.2), 
an aperture of 3 arcsec is equivalent to central regions of ~ 3.5 
kpc. On the other hand, an aperture of the 5 arcsec used to de¬ 
fine the central region in this study covers a central portion of 
~ 1 kpc at the median redshift of our sample (z ~ 0.01). As 
a consequence, the SDSS-pair single-fiber observations cover a 
wider region than the area covered by our central aperture. In 
fact, at the 20 arcsec aperture where the differences in metallic¬ 
ity (using the calibrator from Pettini & Pagel (2004)) are similar 
to those found by Ellison et al. (2008), this aperture covers a 
region of ~ 3.5 kpc at the median redshift of our sample. Al¬ 
though we did not use the same indicators for metallicity as for 
the single-fiber spectroscopic survey studies, we suggest that by 
using abundance calibrators presented in the literature, it is pos¬ 
sible to obtain the same lower metallicity in interacting galaxies. 

Finally we want to highlight that regardless of the indicator 
used to derive the oxygen abundance, the metallicity in the cen¬ 
tral regions for interacting galaxies is similar to the one derived 
for the control sample. 

4. Conclusions 

We used the spatially resolved information provided by the 
CALIFA survey to carry out the first statistical study of the im¬ 
pact of the merger event on the sSFR and the oxygen abundance 
on different galactic scales for a sample over 100 galaxies at dif¬ 
ferent stages of interaction. Moreover, this survey allows us to 
homogeneously determine the same properties in a control sam¬ 
ple of 80 non-interacting galaxies. We separated the galaxies in 
both samples between star-forming and non-star-forming galax¬ 
ies (see Secs. 3.1 and 3.2). For the first subsample, we find a 
(moderate) enhancement in the sSFR in the central region of in¬ 
teracting galaxies. However, in outer regions, the sSFR is sim¬ 
ilar to or moderately suppressed in comparison to the control 
sample (see Sec. 3.2). This agrees with previous observational 
studies (e.g., Ellison et al. 2008; Yuan et al. 2012), as well as 
with numerical simulations (e.g., Hopkins et al. 2013; Moreno 
et al. 2015). These studies indicate that tidally induced star for¬ 
mation becomes evident in the central region of the galaxies as 
a consequence of gas inflows. In this scenario, the new supply 
of metal-poor gas to the central region produces a dilute metal¬ 
licity. In contrast to this picture, we find similar metallicities in 
the central region of star-forming, interacting galaxies to those 
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Fig. 8. Differences in the median values of AO/H between the interacting and control samples at different aperture sizes. Left: The aperture sizes 
are measured in arcsecs. Right: Aperture scales are measured in effective radius. The dashed line represents the difference reported by Ellison 
et al. (2008). Errors bars are determined using bootstrapping in both samples. On any of the scales, the differences reach the values found in 
single-fiber spectrocopic studies. 


X 

0 

T3 


0 

O 

C 

0 

0 


0.04 
0.02 
0.00 
- 0.02 
-0.04 
-0.06 

5 10 15 20 25 30 

aperture diameter (arcsec) 




• Marino et. al (2013) 

♦ Pettini & Pagel (2004) 




Fig. 9. Similar to Fig. 8. Differences in the median values of AO/H 
between the interacting and control samples at different aperture sizes. 
In this case we compare using apertures in arcsec units with the differ¬ 
ences in metallicity using different calibrators. Blue circles represent 
the differences obtained from the 03N2 revisited calibrator used in this 
study given in Marino et al. (2013). Red diamonds show the difference 
obtained from the calibrator presented by Pettini & Pagel (2004). The 
dashed line represent the difference reported by Ellison et al. (2008). 


derived for isolated galaxies (see Sec. 3.3). When we consider 
larger apertures (in effective radius units), we find hints of di¬ 
lute metallicities in the interacting galaxies (see Sec. 3.4). Our 
results support the notion of a tight interplay between different 
physical processes in the central part of interacting galaxies. Al¬ 
though metal-poor gas inflows can be considered as the main 
process that affects the chemical evolution toward the center of 
interacting galaxies, there are some other processes to be taken 


into account that could enrich the ionized gas, such as stellar (or 
AGN) feedback or returned material into the ISM. The results 
presented here encourage IFS studies in larger samples of inter¬ 
acting galaxies in order to understand the evolution of the galac¬ 
tic chemical content and to quantify how different merging con¬ 
figurations can affect the merger-induced metallicity dilution. 
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